Spatial two-particle NOON-states in periodically shaken three-well potentials 
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Few-particle dynamics in a three-well potential are investigated numerically. It is shown that 
periodically shaking the potential can considerably increase the fidelity of emerging spatial quantum 
superpositions. Such NOON-states are important for quantum interferometry. If the two particles 
initially sit in the middle well, the probability to return to this state can distinguish pure quantum 
dynamics from statistical mixtures. The numeric implementation of decoherence via particle losses 
shows clear differences from the pure quantum behaviour. A three-well lattice could be an ideal 
system for experimental realisations. 
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I. INTRODUCTION 

Double- well lattices [UH] provide an interesting, exper- 
imentally accessible tool to investigate quantum effects 
for atom numbers even down to less than six atoms [3J. 
Three- well lattices and four-well lattices could be realised 
either by using subwavelength lattices [4] or by adjusting 
the technique used in Refs. [3J [5J by using more lasers. 
The number of particles in one of the two wells of a 
double-well lattice can easily be measured by averaging 
over all double wells [3J . For a single three- well potential, 
transport and interaction blockade of cold bosonic atoms 
have been discussed in Ref. jl| . In a single three- well po- 
tential, few particle dynamics could be observed by using 
the single-site addressability of Refs. [8] . 

One of the fascinating aspects of quantum mechan- 
ics are non-classical quantum states like the NOON- 
states [PJ: 
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where \N — n,n) refers, e.g., to N — n particles being 
on the left and n particles on the right of some bar- 
rier but NOON-states can also be investigated in phase- 
space. These states have also been called "Schrodinger 
cat"-states, e.g., in Refs. [ID! ITT]. Like the spin-squeezed 
states of Refs. [T2][T3], such states are relevant to improve 
interferometric measurements |14j . Suggestions on how 
such interesting many-particle quantum superpositions 
might be obtained can be found, e.g., in Refs. [T5TE5] 
and references therein. 

The focus of the present paper lies on the numeric cre- 
ation of spatial two-particle NOON-states in three-well 
potentials with a focus on experimental signatures to ver- 
ify that indeed NOON-states rather than statistical mix- 
tures have been created. To achieve the NOON-states 
via controlled quantum dynamics, we suggest to peri- 
odically shake the three-well potential. This approach 
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to enhance tunnelling in a controlled way is based on 
photon-assisted tunnelling [26H5P1 which was realised ex- 
perimentally for Bose-Einstein condensates in optical lat- 
tices [31] (for the interpretation of the experimental data 
see also [32]). Recent research in periodically driven sys- 
tems include nonlinear Landau-Zener processes [33] and 
transport of bound pairs [3JJ [35] . 

The paper is organised as follows: in Sec. [TTJ we in- 
troduce the Hamiltonian used to model the periodically 
shaken three-well potential as well as our main aim, the 
two-particle NOON-states. Section |III| shows that for 
some parameters, the quantum dynamics of the period- 
ically shaken system can be understood within a sim- 
plified model originally developed for double wells (30] 
which we adopt for the present situation. In Sec. |IV| 
we introduce the two-particle NOON-states and explain 
how they can be distinguished from statistical mixtures. 
Section [V] shows that periodic shaking can considerably 
increase the fidelity, i.e., the quality of the quantum su- 
perpositions. In order to show that the quantum su- 
perpositions can indeed be distinguished from statistical 
mixtures, last but not least decoherence effects are in- 
cluded. 



II. MODEL SYSTEM 

For the description of ultracold atoms in a three-well 
potential, the three-site version of a many-particle Hamil- 
tonian [36] originally developed in nuclear physics [37] is 
used: 
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The operator annihilates (creates) a boson in well j; 
J is the hopping matrix element, 2ft/io is the tilt between 
well 1 and well 3 and 2/i/ii is the driving amplitude. The 
interaction between a pair of particles in the same well is 
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FIG. 1. The three-well potential with initially two atoms in 
the middle well. This situation is an experimentally realistic 
initial condition 42. Experimental investigations would, e.g., 
be possible in three-well lattices [4] (cf. the double-well lattice 
potential with less than six atoms per lattice site of Ref. [3]). 
Periodic shaking is currently established as an experimental 
tool to control the quantum dynamics 31, 32, 43, 44 . 



denoted by U. Quantum dynamics beyond such models 
has been investigated, e.g., by Refs. |51 155PD] [4"T] . A 
systematic sketch of the three-well potential is shown in 
Fig. [I] 

A Fock-basis is useful to describe N atoms in the three- 
well potential: 

\lpN-n-m,n,m) = \N - Tl - m, 71, m) 

= \n,m) (3) 

with m particles in the right well, n particles in the mid- 
dle well and N — m — n particles in the left well. 

The aim of this paper is not to produce a single-particle 
quantum superposition but a two-particle NOON-state: 
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As will be shown in this paper, periodic shaking helps 
on the way to achieve this aim. In order to characterise 
the quality with which this aim is achieved for a given 
normalised wave-function \ip(t)) ((ip(t)\ip(t)) = 1), the 
fidelity, 



(5) 



or more precise the fidelity maximised over all possible 
angles 8 is used: 



fidelity = max (F s , < 8 < 2tt) . 



(6) 



In order to be a useful flag to indicate a NOON-state, 
the fidelity has to be larger than 50%; below this value, 
the state might not be a superposition at all as, e.g., the 
state |2, 0, 0) would also result in a fidelity of 0.5. 

Figure [2] displays an example for which the fidelity is 
well above 50%. The data is plotted as a function of the 




FIG. 2. Fidelity Q as a function of dimensionless time r 
[Eq. Q] for the experimentally realistic initial condition |42| 
of two particles sitting in the middle well of a three- well poten- 
tial. The parameters were chosen such that shaking frequency, 
interaction and tilt are equal: Tiuj = U — h/io = 6 J which cor- 
responds to integer photon-assisted tunnelling [271 129| . The 
shaking amplitude was chosen to be h[i\ = 25. 2J. It has been 
optimised with the simplified model (17 1. Solid grey line: the 



exact numerical solution of the time-dependent model Hamil- 
tonian Q. Dashed black line: the simplified model which can 
be solved via exact diagonalisation. While the fidelity is well 
above 50%, it still remains below 90% (the exact numerics 
reaches the maximum 0.8852 at r = 4.787). 



dimensionless time, 
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Figure[3]shows that for the timescales of Fig.[2]and a large 
parameter range, even without shaking, the fidelity can 
reach values well above 50% (Fig. [3]) but it stays again be- 
low 90% (Fig.[3|. Thus, so far there does not seem to be 
an advantage of using periodic shaking. However, before 
showing that periodic shaking can lead to fidelities well 
above 95%, we show that it is possible to quantitatively 
understand the behaviour displayed in Fig. [2] 



III. SIMPLIFIED MODEL 

In order to understand the behaviour under periodic 
shaking, a simplified version of the Hamiltonian ^ 
can be used. This approach introduced below will be 
useful for large shaking frequencies and not too long 
timescales [3D]. 

Extending the simplified model of Ref. [30] to a three- 
well potential starts with writing the Schrodinger equa- 
tion corresponding to the Hamiltonian ^ in the Fock 
basis Using the short notation of Eq. ([3]), we then 
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with Hq 
states is 



= limj_^o H. The tunnelling between Fock 
then included in a system of first-order 



one-dimensional differential equations for the 



(cf. 



Rcfs. [28, 30 ). The time-dependence enters entirely via 
the exponents introduced in the ansatz it includes 
terms linear in time (either proportional to the inter- 
action or to the tilt) as well as a term proportional to 
the driving amplitude for which the time-dependence is 
a cos(wi). This last term can be transferred into a sum 
of exponentials with linear time-dependence via |45j 



iz cos(u;t) 



oo 



Jk{z)i 



k ^ikujt 



(9) 



FIG. 3. Fidelity in a two-dimensional projection as a function 
of both interaction U and tilt hfio without any shaking. To 
obtain the plotted data, the model Q was solved numerically 
for two particles initially in the middle well and times r with 
< t < 50. for each set of parameters, the maximum is 
plotted. Similar to the periodically shaken case displayed in 
Fig. [2| the fidelity stays below 90% (0.892 at U = 1.3 J, hp = 
—0.5 J and r ~ 15). Higher fidelities will be achieved via 
periodic shaking (cf. Fig. [4| 



introduce the ansatz motivated by the interaction picture 



where the Jk{z) are Bessel functions. This leads to a 
series of relevant frequencies which for two particles read: 
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The model is still equivalent to the full Schrodinger equa- 
(8) tion corresponding to Eq. it now reads: 
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where 



A 1 {t) = -jV2Y J Jk{^)i k eM*Vkt) , (14) 

k 

A 2 {t) = -jV2j2Jk{^ L )^ k eM^kt) , (15) 

k 

A 3 (i) = -J^J i ,.(^)z fe exp^) . (16) 



In general, for given parameters of tilt, interaction and 
shaking frequencies, there will be only one term with slow 
time-dependence [35] . Following the reasoning of the ro- 
tating wave approximation 47 , this term should domi- 



nate the dynamics for not too large timescales. Defining 
U 



OV = -r- + 2^o - k'ui 

n 



(17) 



where k' is the integer for which | ¥ + 2fiQ — kuj\ reaches 
its minimum and repeating the same reasoning for 
and cx/c leads to a much simpler version of the matrix 
in Eq. (131 which surprisingly can even be solved via 
exact diagonalisation [15] (see Ref. [3D] where even some 
analytical results could be derived). The derivation of 



our simplified model explains [cf. Eq. (17)] why the two 
curves in Fig. [2] agree well except for the deviation for 
larger times. 
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IV. HOW TO DISTINGUISH TWO-PARTICLE 
NOON-STATES FROM STATISTICAL 
MIXTURES 

The main idea is to be able to distinguish quantum su- 
perpositions from statistical mixtures via the behaviour 
concerning the return to the initial state. Why a tilted 
three-well potential is useful for considerations like that 
can already be motivated on the single-particle level: On 
the one hand, despite the tilt the tunnelling processes 

|0, 1,0) ->-^(|l,0,0)±|0,0,l))->|0,l,0> (18) 

do not involve any change of the potential energy. On 
the other hand, the individual tunnelling processes 



right well at time t\ < t% - really is unambiguous. Using 
the time-evolution operator U, one has: 



|0,1,0) 



11,0,0) 
|0,0,1) 



|0,1,0). 



(19) 



are related to a change in energy (because of the tilt) 
and thus both the first and the second tunnelling pro- 
cess of Eq. (19) would thus be suppressed. This of- 



fers a starting point to distinguish the time-evolution 
of quantum superpositions from statistical mixtures: On 
the one hand, tunnelling of the quantum superposition 
^ (| 1,0,0) ± |0,0, 1)) back to the initial state |0, 1,0) is 
allowed energetically. On the other hand, for a statistical 
mixture like 

q= i(|l,0,0)(l,0,0| + |0,0,l)(0,0,l|) (20) 

the tunnelling to the state |0, 1,0) is suppressed for en- 
ergetic reasons. 

There is, however, an even stronger reason why the 
return to the initial state can verify interesting quantum 
superpositions: 

While experimentally measuring the fidelity ^ would 
either be impossible or at least be a considerable chal- 
lenge for future experiments, measuring the number of 
atoms in one of the wells is comparatively straight- 
forward both in an optical superlattice |3j and in a sin- 
gle three- well potential p3IB]. Thus, for the two-particle 
NOON-state -3= (|2, 0, 0) + exp(iS)\Q, 0, 2)) one could ex- 
perimentally verify, that with nearly 50% probability 
both particles are in the left well and with 50% prob- 
ability both particles are in the right well. In order to 
distinguish this situation from the statistical mixture, 

e=i(|2,0,0)<2,0,0| + |0,0,2)<0,0,2|) , (21) 

we will be searching for parameters such that the proba- 
bility to return to the initial state is high. As the initial 
state is the Fock state with two particles in the middle 
well, this is an unambiguous experimentally verifiable sig- 
nature. 

It remains to be shown that the return to the initial 
state at time t 2 - in combination with the 50/50 prob- 
ability to find both particles either in the left or in the 



m 2 )) = u(t 2 ,t 1 )m 1 )) . 
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Let us consider one particular ideal set of parameters for 
which the wave-function ipi(t) reaches the perfect cat- 
state Q at t = t\ followed by a probability to return to 
the initial state at t = t 2 of 100%: 



|(0,2,0|Vi(t 2 ))| S 



1 



(23) 



Within this ideal set of parameters, an inaccurate prepa- 
ration of the wave-function could lead to a return- 
probability of just x: 

\(Mh)Hi(t 2 ))\ 2 = x , Q<x<l. (24) 

The fact that U^U = 1 implies: 

\(Mti)\Mh))\ 2 = x . (25) 

Thus, for our example of a perfect return to the initial 
state, a fidelity of x at t = t\ leads to a probability of 
x to return to the initial state at t = t 2 . In particular, 
if a NOON-state has a 100% return probability to the 
initial state at a given time, the corresponding statisti- 
cal mixture ( |21[ ) will only manage a return to the initial 
state with 50%probability at the same point of time [15] • 
While we will investigate one particular source of de- 
coherence - particle losses - in the next section, other 
sources of decoherence would also lead to statistical mix- 
tures and thus to a different dynamics. Given the level 
at which today's experiments with ultra-cold atoms are 
performed, we can assume that the particles do mostly 
behave as predicted by the Schrodinger equation without 
decoherence. 

If the return to the initial state is less than the expected 
100% (but still well above 50%), this does not necessar- 
ily imply that the NOON-state was partially destroyed 
via decoherence: experimental uncertainties which would 
average over different solutions of the Schrodinger equa- 
tion could produce a similar effect. Nevertheless, such 
a high return probability can still be a signature of the 
NOON-state: It has to be checked if the variations of the 
parameters still lead to high fidelities at t = t±. In an 
experiment this would lead to probabilities of the order 
of 50/50 to find the two particles either in the left or in 
the right well. 



V. NOON-STATES WITH FIDELITIES 
ABOVE 90% 

To obtain better results than displayed in Figs. [2] and 
[3j three requirements were used for the numeric optimi- 
sation of the parameters: 

• a fidelity well above 90%. 

• a high probability to return to the initial state. 
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FIG. 4. Fidelity Q as a function of time for two particles 
initially sitting in the middle well and hui — 6J, U = 9.1J 
in an untilted three-well potential (no = 0) with a shaking 
amplitude of h[ii — 18. 3J. Black solid line: numeric solution 
of the model j2j; the fidelity reaches values of 99.67%, well 
above 90%. Gray solid line: the simplified model ( |17[ l does not 
even predict the correct frequency of the oscillations; taking 
15 terms in Eqs. ( 14 H 16 1, again leads to a better agreement 
(dashed black line) . 



not too large timescales. 



There are parameter-sets with and without tilt that fulfil 
these conditions. One example is shown in Figs. [I] and [5] 
The approach via the simplified model (17) only works 

in _ 



if more than one frequency is included in Eqs. (|14j)-([16j) 
(see Fig. [4]). At first glance, the standard deviation dis- 
played in Fig.[5jmight render it difficult to distinguish the 
behaviour including decohcrcncc from the ideal quantum 
case. However, for the accuracy with which the average 
can be measured, the standard error is relevant which 
is smaller by a factor of \j\f~M for M measurements. 
Taking M ~ 100 this would thus reduce the error by a 
factor of 10. Our theory so far ignores decoherence. For 
the present system, atom losses via scattering with back- 
ground atoms are an important decoherence mechanism 
(cf. Refs. ED]). To include this into our model, we use 
a method called piecewise deterministic processes, PDP 
in Ref . [SI] , which is based on Refs [S3 [S3] . The main 
idea is to use quantum dynamics for the time-evolution 
in combination with atom losses for which the points 
of time of the atom losses are chosen via random num- 
bers [SB [S3]. The results displayed in the middle panel of 
Fig. [5j are averaged over many sets of random numbers. 
Atoms are lost with a loss-rate a (cf. [A| , 
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FIG. 5. Average number of particles in the second well as a 
function of time for the same parameters as in Fig. [4] Up- 
per panel: thick solid line: (712) as predicted by quantum 
dynamics obtained by numerically solving the model at 
r ~ 13.36 there are hardly any particles in the middle well 
((712) = 0.026) and at r ~ 26.69 the particles have returned 
to the middle well ({TI2) = 1.9919). Dashed lines: (712) ± A712 
where A712 is the standard deviation. Middle panel: thick 
solid line: (712) as predicted by a combination of quantum 
dynamics and particle losses ( |26[ ) for a = J/ (30ft). Black 
dashed lines: (712) ± Ari2 where A712 is the standard devia- 
tion. Lower panel: (712) as predicted by quantum dynamics 
obtained by numerically solving the model |2| if the quantum 
superposition turns into a statistical mixture at the point of 
highest fidelity (r max = 13.356). The required accuracy of the 
parameters is discussed above Eq. ( 28 1 . 



corresponding to an exponential decay of the average 
number of particles, which for the two particles reads: 



(N(t)) = 2exp(-crf) . 



(27) 
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(N(t)) = -a(N(t)} 



(26) 



The middle panel of Fig. [FJ uses a = <X\ with a\ — 
J/(30h); for the timescales of the double- well lattice 
experiment [3 , the fidelity reaches its maximum at 
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mean[|<n 2 >-<n 2 > jd | ] 



with TLp+.ti as the characteristic timescale: 
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FIG. 6. Mean square deviation (29 1 for T = 100ft/ J in two 



dimensional pr ojec tion as a function of both settling time 
Tsctti [see Eq. ( 28 1] and initial phase. All other parameters 
as in Fie [4] T dri vc = 2tt/ ui is the oscillation period. Contrary 
to experiments like [55] , in our case the initial phase does not 
play any role. 



9.7... 16.1ms and thus a « 1/36 ... 1/22 kHz. In or- 
der to be able to observe the return to the initial state 
with 1% accuracy, we require the average number of par- 
ticles (27) to be at least 99% of the initial value up to 



a timescale for which pure quantum evolution predicts a 
nearly perfect return to the initial state (32 ms). This 
results in a lifetime of the particles of l/a^ ~ 3 s. This 
realistic lifetime of a single particle in the lattice lies two 
orders of magnitude below the lifetime required to pro- 
duce NOON-states of, say, 100 particles according to the 
scheme of Ref . [22]. 

The lower panel of Fig. [5] shows what would happen 
if the quantum superposition was turned into a statis- 
tical mixture at the time for which the fidelity reaches 
its maximum: although now there still are two particles 
involved in the simulation, the average occupancy of the 
middle well would remain a factor of two below the values 
reached for pure quantum evolution. 

In order to reproduce the generation of NOON-states 
shown in Figs. [4] and [5j an experiment would require a 
control of both interaction (£// J) and driving amplitude 
hfii of the order of several per cent, the tilt 2fi,/io should 
be zero to about O.QAJ/h in order to reach fidelities above 
0.95 [55]. In an experiment, there might be problems in 
realising a periodic shaking which is exactly zero for times 
t < and perfectly sinusoidal for times t > 0. To show 
that this is not necessarily a problem, Fig. [6] shows what 
happens if the shaking amplitude is switched on and off 
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x 1 — exp {- 
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(28) 

the shaking now includes a phase <f> in sm(ut+<j)). Plotted 
is the mean deviation from the case without any increase 
of the amplitude: 



mean 
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where (n 2 )id is the model 



(n 2 )id\ dt 

(29) 



VI. CONCLUSION 

To conclude, we have shown via numerical investi- 
gations that high-fidelity-two-particle NOON-states can 
emerge in three-well potentials. There is, in particular, 
an easily observable difference between statistical mix- 
tures and perfect NOON-states: only the latter would 
return to the initial states. For realistic conditions, de- 
coherence would be present in such experiments, which 
would destroy the quantum superpositions. Here, the 
proposed scheme might offer the possibility to measure 
the timescales of decoherence. 

Three-well lattices could be used to realise such ex- 
periments - produced either via a superposition of lasers 
similar to the double-well lattices of Refs. [5] [5J or by 
using subwavelength lattices [3]. The number of atoms 
per three-well lattice can be controlled precisely via the 
Mott- insulator used to load the lattice [3]. Such a setup 
would allow to measure the average number of particles 
of one of the three lattice sites in a single experiment 
by averaging over all three-well potentials. Furthermore, 
by focusing on one single three-well potential within the 
three- well lattice via techniques similar to Refs. [7] [5] and 
repeating the experiments many times, even more precise 
measurements could be possible. 
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Appendix A: Lindblad equation to model particle 
losses 



A derivation of the Lindblad Master equation can be 
found, e.g., in Ref. [51 . Often, atf> is the annihilation 
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(creation) operator of a photon in mode u> but atom losses 
can be described by the same approach: 



Pnif) 



(A2) 



K + l t ■■ i - t 
— — (aa 1 g + gaa' 



2apa^) 



2a' ga) 



(Al) 



where p is the density matrix. We assume that particles 
only are lost, thus k + — and a = k_ and thus derive 
a rate equation [?7] by using the fact that the photon 
number distribution p n (t) is related to the density matrix 



The rate-equation thus reads: 
d 



dt 



p n {t) = a(n + l)p n+1 (t) - anp n (t) 



(A3) 



Using this rate-equation, it is straight-forward [47] to 
derive the exponential decay of the numbers of parti- 
cles given in Eq. (26). As explained at the beginning 
of Secjvj_ we use piecewise deterministic processes [5T] 
(cf. [521 153j) to model the atom losses; in between loss 
events, the Hamiltonian dynamics of Eq. ^ is used. 
Further theoretical investigations of atom losses can be 
found, e.g., in Ref. [53] . 
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